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Q ' Abstract 

o : 

^SJ . The main object of our study is a four dimensional Lie algebra which de- 

\ scribes the symmetry properties of a nonlinear Black-Scholes model. This 

■ model implements a feedback effect which is typical for an illiquid market. 

The structure of the Lie algebra depends on one parameter, i.e. we have 
to do with a one-parametric family of algebras. We provide a classification 
of these algebras using Patera-Winternitz method. Optimal systems of 
p one-, two- and three- dimensional subalgebras are described for the family 

, of symmetry algebras of the nonlinear Black-Scholes equation. The opti- 

mal systems give us the possibility to describe a complete set of invariant 
solutions to the equation. 
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<N : 1 Introduction 

^ . In [2] Frey and Patie examined the feedback effect of the option replication strat- 

(3 ! egy of the large trader on the asset price process. They obtained a new model 

by the introduction of a liquidity coefficient which depends on the current stock 
price. The feedback-effect described leads to a nonlinear version of the Black- 
' Scholes partial differential equation, 

a , 

with S G [0, oo), t G [0, T]. As usual, 5* denotes the price of the underlying asset 
and u{S,t) denotes the hedge-cost of the claim with a payoff h{S) which will be 
defined later. The hedge-cost is different from the price of the derivatives product 
in illiquid markets. In the sequel t is the time variable, a defines the volatility of 
the underlying asset. The Lie group analysis of the equation ([T]) was provided in 
[1]. By using this method the Lie point symmetries, the Lie symmetry algebras 
and groups to the corresponding equations were found; see for details [T] and 
[H El [3] for a general introduction to the methodology. 
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Theorem 1.1 (Bordag, [T]). The differential equation (J\) with an arbitrary func- 
tion X{S) possesses a trivial three dimensional Lie algebra Dif f/^{M) spanned by 
infinitesimal generators 

d ^ d d 

ot ou ou 

Only for the special form of the function \{S) = uS^, where uj,k E R equation 
(J\) admits a non-trivial four dimensional Lie algebra L spanned by generators 

with commutator relations 

[vi,V2] = [vi,V3] = [vi,V4] = [V2,V3] = 0, 

[V2, Va] = -kV2, [V3, W4] = (1 - k)v3. (3) 

In the paper ^ authors try unsuccesses to construct an optimal system of 
subalgebras for the symmetry algebra (I2l). The authors used the method sug- 
gested in the series of well known papers by P.Winternitz and Co [3, [8j where all 
three and four-dimensional Lie algebras where classified. 

The investigation in [9J contains some misprints and mistakes which demand 
corrigendum. In that paper the structure of the optimal system of subalgebras 
[9] does not contain some one-dimensional algebras. On the other sides the clas- 
sification does not depend on the parameter k from which the structure of the 
algebra deeply depends. To be able to construct correct families of invariant 
solutions we need a correct optimal system of subalgebras. 

In our paper we present the correct optimal system of one-, two-, three- dimen- 
sional systems of subalgebras. 



2 Classification of the algebra L 

Let us consider the four-dimensional Lie algebra L ([2]) with the commutator re- 
lations To classify this algebra we use the classification method which was 
introduce by J. Patera and P. Winternitz in [7|. Further we use the notations of 
this paper. 

As we noticed before the structure of the algebra L depends on the one real- valued 
parameter k. As it was remarked in |l] the algebra posses a two-dimensional cen- 
ter by A; = and k = 1. 



Case k = 0. In this case the generators of the algebra L take the form 

d d d d d 

Vl = ^, V2 = X—, V3 = TT^ ^4 = X— + U — . (4) 

ot OU Ou Ox ou 
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with the following commutator Table [T] 
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Table 1: The commutator table of the algebra L in case k = 
Let us introduce an algebra L'^ ^ with operators 

ei = -W4, 62 = Vs, 63 = Vl, 64 = V2. (5) 

Then L'^^ = A2(B 1A\ with one nontrivial commutator relation [61, 62] = 62. 



Case k = \. This case leads to the algebra L4.2 with generators 

d d d d 

Vl = TTT, V2 = X — , V3 = ^4 = X — . (6) 

at ou ou ox 

which is isomorphic to the algebra L'^ 2 spanned on generators 

61 = f4, 62 = V2, 63 = Vl, 64 = f3. (7) 

In this case algebra L4 2 has the decomposition A2 © 2Ai with one nontrivial 
commutator relations [61,62] = 62. 



Case > 7^ 1- In this case the algebra L is isomorphic to the algebra 
L43 = A3 5 © Ai with the commutator relations [61,63] = 61, [62,63] = ae2, 
where ^ 

61 = V2, 62 = W3, 63 = --^va, 64 = Vl (8) 

and a = 

k 



Case A; < |, A; ^ 0. This case leads us to the algebra L which is isomorphic 
to L4.4 = yig g © Ai with the commutator relations [61,63] = 61, [62,63] = 062, 
where 

1 

ei = V3, 62 = V2, 63 = YTr^^4, 64 = Vl (9) 

and a = . 
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Case k = \. In the last case the algebra L has generators of the following type 
d d d did 

Vl = V2 = X — , V3 = Vi = X— + -U — , (10) 

ot ou ou ox 2 ou 

and the commutator Table [2] 





Vl 


V2 


vs 




Vl 
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V3 













V4 







-|^3 






Table 2: The commutator table of the algebra L in case k = ^ 

and is isomorphic to the algebra L'^ ^ = A^ ^ © Ai where generators are denoted 
by 

ei = ^^3, 62 = V2, 63 = 2Vi, 64 = Vi- (11) 

The nontrivial commutator relations are [61, 63] = 61, [62, 63] = —62. 



3 An optimal system of subalgebras 

The main goal of this paper is to find a correct optimal system of subalgebras for 
the Lie algebra L ([2]). The procedure was described by Pattera & Winternitz in 
[7j. In the paper all three and four-dimensional algebras were classified and the 
optimal systems for these algebras were listed. We repeat this algorithm for the 
algebra L first in general case where k ^ 0,^,1. In those two cases ([H]), (E]) the 
Lie algebra L is isomorphic to the algebra g © Ai with following commutator 
Table [3] 





ei 


62 


63 


64 


ei 








ei 





62 








a62 





63 


-ei 


-ae2 








64 















Table 3: The commutator table of the algebra L where < |a| < 1. 

In the general case L has one central element and can be represent as a direct 
sum of one- and three- dimensional Lie algebras 

L = {64}©L3, (12) 
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where 64 is the central element of the algebra L and L3 = L \ {64}. 

The representation (fT2l) simplify the procedure of construction of an optimal 

system of subalgebras. 

We start with construction of the corresponding system of subalgebras for both 
algebras in [12] and then complete the study with the investigation of non-splitting 
extensions. We follow the paper [7] and describe a solution of this problem in a 
step-by-step method introduced by the authors. 

Step 1. We find all subalgebras of {64}, it means we have {0} and {64}. 
Step 2. We have to classify all subalgebras of L3 ( fT2l) under conjugation which 
is defined by an interior isomorphism of this algebra. This isomorphism can be 
presented by the adjoint representation. 

Definition 3.1 (Olver, [6]). Let G he a Lie group. For each g & G, group 
conjugation Kg{h) = ghg~^, h & G, determines a diffeomorphism on G. More- 
over, Kg o Kg/ = Kgg/,Ke = 1 ^ , SO Kg detcrmines a global group action of G 
on itself, with each conjugacy map Kg being a group homomorphism: Kg{hh') = 
Kg{h)Kg{h') etc. The differential dKg : TG\h — ^ *TG\Kg{h) is readily seen to 
preserve the right-in variance of vector fields, and hence determines a linear map 
on the Lie algebra ofG, called the adjoint representation: 

Ad g{v) = dKg{v) (13) 

The simplest way to find the adjoint representation is the Lie series 

Ad{e'x.\){ev))w = w — e[v^ w] + "^bi [^-i w]] — . . . (14) 

The adjoint representation table for the algebra L f|T2|) is rather simple and has 
a form given in Table H] 



Ad 


61 


62 


63 


64 


ei 


61 


62 


63 - eei 


64 


62 


61 


62 


63 — aee2 


64 


63 


e^Ci 


e"^e2 


63 


64 


64 


61 


62 


63 


64 



Table 4: The adjoint representation table of the algebra L = A'^^ © Ai with 
{i,j)-th entry indicate Ad{exp{eei))ej element. 

By using the adjoint representation IHM we classify all subalgebras of L3 (|T2ll 
under conjugacy. 
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One-dimensional subalgebras. Firstly we consider one- dimensional subalge- 
bras of the general type 

A = {aei + be2 + ces}, 

where a, b, c are arbitrary constants. If c 7^ then we can use the first and the 
second lines of the Table HI We obtain 

Ad{exp{^ei + (^2))^ = (a - c^)ei + {b - caC)e2 + ces, (15) 

choosing ^ = -, ( = — we prove that A is isomorphic to 63. 
If c = in (fT5|) we have three cases to study. 

If a 7^ 0, b = then A is isomorphic to ei. If a = 0, 6 7^ then A is isomorphic 
to 62- The last case we obtain if a6 7^ then we can use the third line of the 
Table H] and obtain after an action of the adjoint representation 

Ad{exp{^e3))A = ae^ei + 6e°«e2. (16) 

Using the scaling on ae^ and choosing ^ = log | f | we prove that A is isomor- 
phic to the algebra generated by {ei ±62}. 

Collecting all previous results we obtain that the optimal system of one- 
dimensional subalgebras of L3 contains following subalgebras 

{0}, {ei}, {62}, {es}, {ei±e2}. (17) 

Two-dimensional subalgebras. Let us consider now two-dimensional sub- 
algebras of L3. Let i? be one of the one-dimensional subalgebras (IT7|) and 
A = {061 + 662 + 063}. For a subalgebra M = B + Awe demand that [A, B] C M. 
Let B = {ei} then without loss of generality we can represent A in the form 
{a62 + be^}. Let f> 7^ 0, by using of the second line of the adjoint representation 
Table m we prove that A is isomorphic to {63}. If 6 = then A = {62}. In this 
case we obtain two subalgebras 

{61,62}, {61,63}, (18) 

which are non-isomorphic to each other. In the same way we obtain subalgebras 
{61, 62}, {62, 63} in case B = 62. 

Let B = {63} then without loss of generality we can represent A = {aei + 662}. 
Let us check the commutator relations 

[a6i + be2, 63] = a6i + abe2. 

We see that A + B is an algebra just under condition ab = 0. On this way we 
obtain the two subalgebras 

{61,63}, {62,63} (19) 

In case B = {ei ± 62} we choose A = {63}. Then [61 ± 62, 63] = 61 ± 0:62 is not 
an algebra. 
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The optimal system of the one- and two-dimensional subalgebras of L3 con- 
tains following subalgebras 

{0}, {ei}, {62}, {63}, {eiiea}, {61,62}, {61,63}, {63,63}. (20) 

Step 3. We have to find all splitting extensions of the algebra {64}. To do this 
we have to find all subalgebras A''^ of L3 such that 

[64, Na] C Na (21) 

and classify all such subalgebras under Nor 164. 

Definition 3.2 (Ovsyannikov, [5j). Let N be a suhalgebra of the Lie algebra 
L. By the normalizer Nor^N of N in L we mean the maximal subalgebra of L 
containing N in which N is an ideal, i.e. 

NoriN = {y e L:[y,x\e N ^xeN}. (22) 



As soon as 64 is a central element and Notlc^ = L any adjoint representation 
does not affect on 64 and Na is any subalgebra of L3. This step is trivial and we 
obtain the subalgebras of the type {64, S} where S running through all subalge- 
bras (|20|). 

Step 4. We have to find all subalgebras of type 

64 + ^a,ei,A^a| , (23) 

where Na is a subalgebra of L3 such that Nor^Na is not contained in L3, G M 
are not all equal to zero and the generators 64 + ^ ajCj are not conjugate to 

i 

64^. Since 64 is a central element of L all of those conditions are satisfied. Let 
Na running through the list of algebras (17) and let ^4 = 64 + ^Ojej. Let first 

i 

Na = {0}. This case is trivial because 64 is the central element and the procedure 
was described on the Step 2. We obtain four subalgebras 

{aei + 64}, {062 + 64}, {063 + 64}, {a(6i ± 62) + 64}, a^O (24) 

If we scale by ^ all of these subalgebras we see that for the two first subal- 
gebras we can use adjoint representation generated by ei and 62 to reduce these 
subalgebras to simplest one. We obtain from Table H] following two subalgebras 

{61 + 66-^64}, {62 + 66""^ 64}, (25) 
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where b = ^ ^ 0. Choosing e = log |^| in the first case and £ = ^ log |^| in the 
second one, we finally obtain the following list of one dimensional non-splitting 
extensions 

{ei±e4}, {e2±e4}, {63 + ae^}, {ci ± 62 + 064}, 07^0. (26) 

Let us now consider two-dimensional non-splitting extensions. To simplify these 
procedures we use as Na the subalgebras of the list (l20l) . We notice that under 
action of the adjoint representation the general form of 64 + ^ hold. 

i 

Let Na be equal to {ei} then without loss of generality we can represent A 
as {64 + 0262 + 0363}. If 03 is not equal to zero we can use the second line of the 
adjoint representation Table [Hand reduce the algebra A to 64 + 063. We obtain 
the following subalgebra 

{63 + ae4, ei}, a 7^ 0. (27) 

In the case 03 = we rewrite A = {ae^ + 62} and use the third line of the adjoint 
representation Table H] to obtain 

{e""e2 + ae4,e"ei} (28) 

or 

{62 + ae-°^e4,e(^-")^ei},a ^ 0. (29) 

By choosing £ = ^ log \ a\ and scaling the second generator of the algebra above 
by the corresponding constant we obtain the following algebra 

{62 ±64,61}. (30) 

The same procedure in the case Na = 62 lead us to the non-isomorphic subalgebras 

{61 ±64, 62}, {63 + 064,621,07^0. (31) 

Let us consider the case Na = {63}, then we can choose A = {64 -|- oiei -f- 0262}. 
Note that A + Na is an algebra just in case 0102 = 0. Those subalgebras were 
considered in the previous cases. Let Na = {ei ± 62} then {A^^, 64 + aiCi + 0262 + 
0363} is an algebra only if 03 = 0. Without loss of generality we can represent A 
as an algebra generated by {ae^ + 61} then by using the third line of the Table H] 
we see that algebras {A, Na} are isomorphic to the following algebra 

{61 + 64,61 + 062}, (32) 

where a G M. Note that the case a = we consider on the third step. Finally we 
obtain the following subalgebra 

{61 + 64,061 + 62}, (33) 
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where a 7^ 0. 

Now we consider case Na = {ci, 62} here we can represent A as {ae^ + 63} and 
obtain the following three dimensional subalgebra 

{63 + 064,61,62}, (34) 

where a 7^ 0. It is easy to see that the other choices of A''^ do not provide any 
other non similar subalgebras. 

We obtain the following list of the optimal system of subalgebras to the algebra 

{0}, {ei}, {62}, {63}, {61 + 62}, {61,62}, {61,63}, {62,63}, 
{64}, {61,64}, {62,64}, {63,64}, {61 + 62,64}, {61,62,64}, 
{61,63,64}, {62,63,64}, {61 + 64}, {62 ±64}, {63 + 064}, 
{61 + 62 + 064}, {63 + 064,61}, {62 + 64,61}, {61 + 64,62}, 
{63 + 064, 62}, {61 + 64, 061 + 62}, {63 + 064, 61, 62}. 



Finally we obtain the complete optimal system of subalgebras of Lie algebra L 
(see Table |5]). 



Dimension 


Subalgebras 


1 
2 
3 


{62}, {64}, {ei + 062}, {61 + 664}, 
{62 + 664}, {63 + 064}, {61 + 662 + 064} 
{61,62}, {61,64}, {62,64}, {63,64}, {61 + 662,64} 
{62 + 664,61}, {61 + 664,061 + 62}, {63 + 064,61}, {63 + 064,62} 
{61,62,64}, {61,63,64}, {62,63,64}, {61,62,63 + 064}, 



Table 5: The optimal system of subalgebras of the algebra L ([2]) in case k 7^ 0, |, 1, 
were o G M, e = +1. 
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We remark now that in case k = ^ the structure of the algebra L is the same 
as in the case above hence the optimal system of subalgebras is the same. 
For A; = or A; = 1 we obtain the following system of subalgebras by the similar 
procedure 



Dimension 


Subalgebras 


1 


{62}, {escosip + e^smcp}, 




{ei + a(e3 cos Lp + 64 sin Lp) } , 




{e2 + e(e3 cos(/? + 64 sin y?)} 


2 


{ei + a(e3Cos(/7 + e4sin(^),e2}, {63,64}, 




{61 + a(63 cos tp + sin Lp) , 63 sin — 64 cos (p}, 




{62 + e(63 cos tp + 64 sin 63 sin — 64 cos (y?}. 




{62,63 sin c/? — 64 cos 


3 


{61,63,64}, {62,63,64}, 




{61 + a(e3 cos Lp + 64 sin (/?), 63 sin (/? — 64 cos 62}, 



Table 6: The optimal system of subalgebras of the algebra L ([2]) in case A; = or 
A; = 1 where a G M, e = ±1, 0<ip<U. 
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4 Conclusion 

In this chapter we return to our original algebra L ([2]) and introduce the optimal 
system of subalgebras Table [5] and Table [6] in original generators. 



Parameter 


Generators 


Optimal System 
of subalgebras 


k = 


V3 = -§-, V4,=X-^+U-§- . 
^ Ou 1 ^ dx au 


{113}, cos sin (^}, {i)4+a(i)i cos (^+t'2 sini^)}, 
\vz+t{v\ cos ¥'+f2 sin v')), \v\+a(v\ cos (j3+f2 sin u:>),vz\, 
{v\,V2}, {v4,+a{v\ cos if+V2 sin(^),fi sin ip—V2 cos i^}, 
{f 3+e(fi cos if+V2 sin sin ip—V2 cos yp}, {1)3,111 sin ip—V2 cos i/p}, 
{i)4,fi,'f2}, {f3,i'l,^'2}, {i'4+a('fl cosip+i)2 sin^j),?)! sinip— 112 cosipi^'a}- 


k = l 


fl = ^i ^2=x-i-, 

at^ au ' 

'^3 = £, Vi=X^. 


{^2}, {vi cos ip+vs sin {fi} , {u4+a(i)i cosv'+^'a sin (/?)}, 
{v2+e{vi cos 1/3+113 sin (/?)}, {i)i+a(i)i cos ip+V3 sin i/j),i;2}, 
{''I I'l'sli {f4+a(f 1 cos (/p+fa sini/)),?;! sin i/J— 1)3 cos i/p}, 
{d2 +6(1)1 cos <p+V3 sin i/j),vi sin tp—v^ cos 75}, {1)2,1)1 sin 1/5—1)3 cos ip}, 
{v4,,vi,V3}, {vi,V2,V3}, {D4+a(t)i cos (/3+i)3 sin (/;), 1)1 sin 1/3—1)3 cos (/),i;2}. 


k < ^, ky^O 


Vl = -§-,, V2=X^, V3 = ^, 


{1)2}, {di}, {D3+ai'2}, {d3+<:i'i}, {d2+':di}, {D4+ai'l}, 
{D3+eD2+aDi}, {1)3,1)2}, {1)3, Di}, {d2,di}, {d4,di}, {D3+ei'2,Di} 
{D2+eDl,D3}, {D3+eDi,aD3+D2}, {D4+aDl,D3}, {D4+aDl,D2}, 
{1)3,1)2,1)1}, {d3,D4,di}, {112,1)4,1)1}, {d3,D2,D4+(IDi}- 


k = \ 


vi = -§-,, V2=x£, V-i = ^, 

V4=X-0- + ^U-^. 

ax 2 au 


{D2}, {di}, {D3+ai'2}, {d3+<:i'i}, {d2+<:i'i}, {114+01)1}, 
{v3+eV2+avi}, {1)3,1)2}, {d3,i'i}> {1)2,1)1}, {'"4,,vi}, {D3+ei'2,Di} 
{D2+eDi,D3}, {D3+efl ,01)3 +1)2}, {1)4+01)1,1)3}, {D4+aDl,D2}, 
{d3,D2,Di}, {1)3,1)4,1)1}, {1)2,1)4,1)1}, {D3,D2,D4+ai'l}- 


A; > i, k^l 




{1)3}, {di}, {1)2+11)3}, {D2+i;fi}, {D3+<:fi}, {d4+cii'i}, 
{D2+eD3+aDi}, {1)2,1)3}, {d2,di}, {1)3, fi}, {f4,Dl}, {f2+eD3,Di} 
{D3+eDi,D2}, {D2+efi,ai)2+i)3}i {D4+ai)i ,D2}, {D4+aDi,D3}, 
{1)2,1)3, Dl}, {1)2,1)4,1)1}, {d3.1'4,1)i}, {D2,D3i1'4+ai)l}- 



Table 7: The optimal system of subalgebras of the algebra L ([2]) with a G M, 

e = ±1, 0<ip<U 
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